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We study a double quantum dot each dot of which is tunnel-coupled to superconducting leads. In 
the Coulomb blockade regime, a spin-dependent Josephson coupling between two superconductors 
is induced, as well as an antiferromagnetic Heisenberg exchange coupling between the spins on the 
double dot which can be tuned by the superconducting phase difference. We show that the correlated 
spin states — singlet or triplets — on the double dot can be probed via the Josephson current in a 
dc-SQUID setup. 
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In recent years, electronic transport through strongly 
interacting mesoscopic systems has been the focus of 
many investigations In particular, a single quantum 
dot coupled via tunnel junctions to two non-interacting 
leads has provided a prototype model to study Coulomb 
blockade effects and resonant tunneling in such systems. 
These studies have been extended to an Anderson im- 
purity 1^ or a quantum dot coupled to superconductors 
[3-^. In a number of experimental [|j and theoretical 
[i| papers, the spectroscopic properties of a quantum dot 
coupled to two superconductors have been studied. Fur- 
ther, an effective dc Josephson effect through strongly 
interacting regions between superconducting leads has 
been analyzed |^|-^. More recently, on the other hand, 
research on the possibility to control and detect the spin 
of electrons through their charges has started. In par- 
ticular in semiconducting nanostructures, it was found 
that the direct coupling of two quantum dots by a tun- 
nel junction can be used to create entanglement between 
spins , and that such spin correlations can be observed 




FIG. 1. Upper panel: sketch of the superconductor-double 
quantum dot-superconductor (S-DD-S) nanostructure. Lower 
panel: schematic representation of the quasiparticle energy 
spectrum in the superconductors and the single-electron levels 
of the two quantum dots. 



in charge transport experiments 1 11 1 . 

Motivated by these studies we propose in the present 
work a new scenario for inducing and detecting spin cor- 
relations, viz., coupling a double quantum dot (DD) to 
superconducting leads by tunnel junctions as shown in 
Fig. It turns out that this connection via a supercon- 
ductor induces a Heisenberg exchange coupling between 
the two spins on the DD. Moreover, if the DD is arranged 
between two superconductors (see Fig. H), we obtain a 
Josephson junction (S-DD-S). The resulting Josephson 
current depends on the spin state of the DD and can be 
used to probe the spin correlations on the DD. 

Model — The double-dot (DD) system we propose is 
sketched in Fig. |^: Two quantum dots (a, 6), each of 
which contains one (excess) electron and is connected to 
two superconducting leads (L, R) by tunnel junctions (in- 
dicated by dashed lines) Q . There is no direct coupling 
between the two dots. The Hamiltonian describing this 
system consists of three parts, Hs + Hd + Ht = Hq + Ht- 
The leads are assumed to be conventional singlet super- 
conductors that are described by the BCS Hamiltonian 

= E / ^ I E i^Ur)h{r)Mr) (1) 
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where ilj is the volume of lead j, h{r) = {—ifiV 



f A)^/2m - and Aj(r) = A^- 



i0j(r) 



is the pair po- 



tential. For simplicity, we assume identical leads with 
same chemical potential /i, and A^ = A/j — A. The two 
quantum dots are modelled as two localized levels and 
£(, with strong on-site Coulomb repulsion J7, described by 
the Hamiltonian 
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where we put = Cfc = ^ e (e > 0) for simplicity. U is 
typically given by the charging energy of the dots, and 
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we have assumed that the level spacing of the dots is '--^ U 
(which is the case for small GaAs dots 0), so that we 
need to retain only one energy level in Ho- Finally, the 
DD is coupled in parallel (see Fig. |l|) to the supercon- 
ducting leads, described by the tunneling Hamiltonian 
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where rj^n is the point on the lead j closest to the dot 
n. Here, $0 = hc/2e is the superconducting flux quan- 
tum. Unless mentioned otherwise, it will be assumed that 
rL,a = rL,b = rL and r_R,a = r_R,h = r_R,. 

Since the low-energy states of the whole system are 
well separated by the superconducting gap A as well as 
the strong Coulomb repulsion U (A, e ^ U — e), it is 
sufficient to consider an effective Hamiltonian on the re- 
duced Hilbert space consisting of singly occupied levels 
of the dots and the BCS ground states on the leads. To 
lowest order in Ht, the effective Hamiltonian is 



H,ff ^PHt [{Eo - Ha)-\l - P)HTf P 
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where P is the projection operator onto the subspace 
and Eq is the ground-state energy of the unperturbed 
Hamiltonian Hq. (The 2nd order contribution leads to 
an irrelevant constant.) The lowest-order expansion (^) 
is valid in the Hmit F < A, e where F = 'Kt^N{Q) and 
iV(0) is the normal-state density of states per spin of the 




FIG. 2. Partial listing of virtual tunneling processes con- 
tributing to H^ff (^. The numbered arrows indicate the 
direction and the order of occurrence of the charge transfers. 
Processes of type (a) and (b) give a contribution proportional 
to Jo, whereas those of type (c) and (d) give contributions 
proportional to J. Other processes not listed here give negli- 
gible contributions in the energy regions of interest. 



leads at the Fermi energy. Thus, we assume that F <^ 
A, e <C U — e, and temperatures which are less than e 
(but larger than the Kondo temperature). 

Effective Hamiltonian — There are a number of vir- 
tual hopping processes that contribute to the effective 
Hamiltonian (|j), see Fig. || for a partial listing of them. 
Collecting these various processes, one can get the effec- 
tive Hamiltonian in terms of the gauge-invariant phase 
differences and Lp between the superconducting leads 
and the spin operators Sa and S;, of the dots (up to a 
constant and with h = I) 



Heff = Jo cos(7r/AB) cos{(l) - irfAB) 

+ [(2 Jo -f J)(l -I- cos(/3) + 2Ji(I + cost: Jab)] 
X [S, • Sfc - 1/4] . 



(5) 



Here Jab = ^ab/^o and ^ab is the Aharonov-Bohm 
(AB) flux threading through the closed loop indicated 
by the dashed lines in Fig. |l|. One should be careful to 
define gauge-invariant phase differences and (/? in (H). 
The phase difference is defined as usual |13| by 



<t> = 0l(i"l) - (t)R{rR) ■ 
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where the integration from to r^^ runs via dot a (see 
Fig. 1^). The second phase difference, ip, is defined by 



= 4>l{vl) - (l)R{rR) 
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The distinction between (j) and cp, however, is not signif- 
icant unless one is interested in the effects of an AB flux 
through the closed loop in Fig. |I] (see Ref. |]ll| for an ex- 
ample of such effects) . The coupling constants appearing 
in (ph are defined by 
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g{xYg{y)[g{x)+g{y)][f{x) + f{y)] ' 



where ( = e/A, f{x) = vT+a?, and g{x) = Vl -I- x^+C- 
Eq. (^) is one of our main results. A remarkable fea- 
ture of it is that a Heisenberg exchange coupling between 
the spin on dot a and on dot b is induced by the super- 
conductor. This coupling is antiferromagnetic (all J's are 
positive) and thus favors a singlet ground state of spin a 
and b. This in turn is a direct consequence of the assumed 
singlet nature of the Cooper pairs in the superconductor 
[Q. As discussed below, an immediate observable con- 
sequence of Heff is a spin- dependent Josephson current 
from the left to right superconducting lead (see Fig. |l|) 
which probes the correlated spin state on the DD. 

The various terms in (^) have different magnitudes. In 
particular, the processes leading to the Ji term involve 
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quasiparticles only as can be seen from its AB-flux depen- 
dence which has period 2$o- In the hmits we will consider 
below, this Ji term is small and can be neglected. 

In the limit ^ 1, the main contributions come from 
processes of the type depicted in Fig. |^ (a) and (b) , mak- 
ing Jo ~ 0.1(r^/Ce) logC dominant over J and Ji. Thus, 
can be reduced to 



Jo cos{tt /ab) cos(0 - TrfAs) 



Heff ~ Jo 

+ 2Jo(l + cos If ) 
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up to order (logC)/C- As can be seen in Fig. g (a), the 
first term in (||) has the same origin as that in the single- 
dot case : Each dot separately constitutes an effective 
Josephson junction with coupling energy — Jo/2 (i.e. tt- 
junction) between the two superconductors. The two re- 
sulting junctions form a dc SQUID, leading to the total 
Josephson coupling in the first term of (^. The Joseph- 
son coupling in the second term in (|^), corresponding to 
processes of type Fig. || (b), depends on the correlated 
spin states on the double dot: For the singlet state, it 
gives an ordinary Josephson junction with coupling 2 Jq 
and competes with the first term, whereas it vanishes for 
the triplet states. Although the limit A^e<C?7 — eis 
not easy to achieve with present-day technology, such a 
regime is relevant, say, for two atomic impurities embed- 
ded between the grains of a granular superconductor. 

More interesting and experimentally feasible is the case 
C ^ 1. In this regime, the effective Hamiltonian ^ is 
dominated by a single term (up to terms of order 



J(l + COS (f) 



(10) 



with J w 2rVe- The processes of type Fig. || (b) and (c) 
give rise to (|0|). Below we will propose an experimental 
setup based on (p^). 

Before proceeding, we digress briefly on the depen- 
dence of J on the contact points. Unlike the processes 
of type Fig. || (a), those of types Fig. ^ (b), (c), and (d) 
depend on 5rL = \rL,a-YL.b\ and Jr^ = |rfl,a - r_R,6|, see 
the remark below Eq. (|^). For the tunneling Hamiltonian 
(^, one gets (putting Sr = Stl — Sru) 



JiSr) ^ ^ 
e 



duj F^{6r,Lj)- F^{6r,uj) 
2tt w + e 
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where F^/^{r,uj) is the Fourier transform of the 
Green's function in the superconductors, F^/^{Y,t) = 
Tje(±t)({V'T(r,i),i/'i(0,0)}) ||. E.g., in the limit e < 
A < /X, we find J{5r) « J(0)e-2«'-/« sin^ikpSr) / {kpSr^ 
up to order l/kp^, with kp the Fermi wave vector in 
the leads. Thus, to have J(Sr) non-zero, 5r should not 
exceed the superconducting coherence length ^. 

Probing spins with a dc-SQUID — We now propose 
a possible experimental setup to probe the correlations 
(entanglement) of the spins on the dots, based on the 




FIG. 3. dc-SQUID-like geometry consisting of the S-DD-S 
structure (filled dots at the top) connected in parallel with 
another ordinary Josephson junction (cross at the bottom). 

effective model (|lo|). According to ([l0| ) the S-DD-S struc- 
ture can be regarded as a spin- dependent Josephson junc- 
tion. Moreover, this structure can be connected with an 
ordinary Josephson junction to form a dc-SQUID-like ge- 
ometry, see Fig. ^. The Hamiltonian of the entire system 
is then given by 



H = J[l + COS(0 - 27Tf)] (Sa • Sb - i 

-I- aJ(l - cos6) , 



(12) 



where / = ^/^q, $ is the flux threading the SQUID 
loop, 6 is the gauge-invariant phase difference across the 
auxiliary junction (J'), and a = J' / J with J' being the 
Josephson coupling energy of the auxiliary junction [ p^ . 
One immediate consequence of ( |T^ ) is that at zero tem- 
perature, we can effectively turn on and off the spin ex- 
change interaction: For half- integer flux (/ = 1/2), sin- 
glet and triplet states are degenerate aX 9 = 0. Even at 
finite temperatures, where 6 is subject to thermal fluc- 
tuations, singlet and triplet states are almost degener- 
ate around 9 = 0. On the other hand, for integer flux 
(/ = 0), the energy of the singlet is lower by J than that 
of the triplets. 

This observation allows us to probe directly the spin 
state on the double dot via a Josephson current across 
the dc-SQUID-like structure in Fig. |[ The super- 
current through the SQUID-ring is defined as Is = 
{27rc/^Q)d{H)/d9, where the brackets refer to a spin ex- 
pectation value on the DD. Thus, depending on the spin 
state on the DD we find 



Is/Ij 
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(singlet) 
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where Ij = 2eJ /Ti. When the system is biased by a dc 
current / larger than the spin- and flux-dependent critical 
current, given by maxe{|/s|}, a finite voltage V appears. 
Then one possible experimental procedure might be as 
follows (see Fig. ^). Apply a dc bias current such that 
alj < I < [a -\- 1)1 J. Here, alj is the critical current 
of the triplet states, and (a -|- l)7j the critical current of 
the singlet state at / = 0, see (|l^). Initially prepare the 
system in an equal mixture of singlet and triplet states by 
tuning the flux around / = 1/2. (With electron g-factors 
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g ~ 0.5-20 the Zeeman splitting on the dots is usually 
small compared with ksT and can thus be ignored.) The 
dc voltage measured in this mixture will be given by (Vo + 
3Vi)/4, where Vb(Vi) ^ 2A/e is the (current-dependent) 
voltage drop associated with the singlet (triplet) states. 
At a later time t ~ 0, the flux is switched off (i.e. / — 
0), with / being kept fixed. The ensuing time evolution 
of the system is characterized by three time scales: the 
time Tcoh ~ max{l/A,l/r} ^ l/F it takes to establish 
coherence in the S-DD-S junction, the spin relaxation 
time Tspin on the dot, and the switching time r^^ to reach 
f — 0. We will assume Tcoh Tspin,Tswi which is not 
unrealistic in view of measured spin decoherence times in 



GaAs exceeding 100 ns Wn. If Ts- 



< T, 



spin 1 



the voltage is 



given by 3Vi/4 for times less than Tspin, i.e. the singlet 
no longer contributes to the voltage. For t > Tspin the 
spins have relaxed to their ground (singlet) state, and 
the voltage vanishes. One therefore expects steps in the 
voltage versus time (solid curve in Fig. Q). If Tgpin < t^^, 
a broad transition region of the voltage from the initial 
value to will occur (dashed line in Fig. m |jl8[. 
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FIG. 4. Schematic representation of dc voltage V vs. time 

when probing the spin correlations of the DD. The flux 

through the SQUID loop is switched from / = 1/2 to / = 



at t = 0. Solid line: Tsw < Ts 



Dashed line: 



To our knowledge, there are no experimental reports 
on quantum dots coupled to superconductors. However, 
hybrid systems consisting of superconductors (e.g., Al or 
Nb) and 2DES (InAs and GaAs) have been investigated 
by a number of groups |19|. Taking the parameters of 
those materials, a rough estimate leads to a coupling en- 
ergy J in (|l^) or (^ of about J - 0.05-0. 5K. This 
corresponds to a critical current scale of Ij ~ 5-50nA. 

In conclusion, we have investigated double quantum 
dots each dot of which is coupled to two superconductors. 
We have found that in the Coulomb blockade regime the 
Josephson current from one superconducting lead to the 
other is different for singlet or triplet states on the double 
dot. This leads to the possibility to probe the spin states 
of the dot electrons by measuring a Josephson current. 
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